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ABSTRACT 

The X-ray pulsar GX 1+4 was observed with the RXTE satellite for a total of 51 ks 
between 1996 July 19 - 21. During this period the flux decreased smoothly from an 
initial mean level of ss 6 x 10 36 ergs -1 to a minimum of « 4 x 10 35 ergs -1 (2-60 keV, 
assuming a source distance of 10 kpc) before partially recovering towards the initial 
level at the end of the observation. 

BATSE pulse timing measurements indicate that a torque reversal took place 
approximately 10 d after this observation. Both the mean pulse profile and the pho- 
ton spectrum varied significantly. The observed variation in the source may provide 
important clues as to the mechanism of torque reversals. 

The single best-fitting spectral model was based on a component originating from 
thermal photons with kTo » 1 keV Comptonised by a plasma of temperature kT » 
7 keV. Both the flux modulation with phase during the brightest interval and the 
evolution of the mean spectra over the course of the observation are consistent with 
variations in this model component; with, in addition, a doubling of the column density 
tin contributing to the mean spectral change. 

A strong flare of duration ^ 50 s was observed during the interval of minimum 
flux, with the peak flux w 20 times the mean level. Although beaming effects are 
likely to mask the true variation in M thought to give rise to the flare, the timing of a 
modest increase in flux prior to the flare is consistent with dual episodes of accretion 
resulting from successive orbits of a locally dense patch of matter in the accretion disc. 
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1 INTRODUCTION 



The study of X-ray pulsars has been an area of active re- 
search for almost 30 years. In spite of this there remain some 
significant shortfalls in the understanding of these objects. 
An example is t he persistent pulsar GX 1+4. At th e time 
of its discovery (Lewin, Ricker and McClintock 1971) it was 
one of the brightest objects in the X-ray sky. The compan- 



ion to the neutro n star is an M6 giant (Davidsen, Malina 



and Bowyer 1977). GX 1+4 is the only known X-ray pulsar 



in a symbiotic system. Measurements of the average spin-up 
rate during the 1970s gave the largest value recorded for any 
pulsar (or in fact any astronomical object) at ~ 2 per cent 
per year. Inexplicably, the average spin-up trend reversed 
around 1983, switching to spin-down at approximately the 
same rate. Since that reversal a number of changes in the 
sign of the torque (as inferred from the ra te of change of 
the pulsar spin period) have been observed ( Chakrabarty et 



al. 1997 ). Several estimates (Beurle et al. 1984, Dotani et 
al. 1989, Greenhill et al. 1993, Cui 1997) indicate a neutron 
star surface magnetic field strength of 2 — 3 x 10 13 G. 

The X-ray flux from the source is extremely variable on 
time-scales of seconds to decades. Two principal flux states 
have been observed, a 'high' state which persisted during 
the spin-up period of the 1970s, and a 'low' state since. Al- 
though the mean flux has been increasing steadily during 
the current 'low' state it has not yet returned to the level 
of the 1970s. Superimposed on these long-term variations 
are smooth changes in the flux on time-scales of order hours 
to days. On the shortest time-scales the periodic variation 
due to the neutron star's rotation period at around 2 min is 
observed. 

Compared to other accretion-powered X-ray pulsars, 
GX 1+4 has a n atypically hard spectrum exten ding out well 
past 100 keV (Frontera and Dal Fiume 1989). Historically 
the spectrum has been fitted with thermal bremsstrahlung 
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or power law models; more recent observations with im- 
proved spectral resolution generally favour a power law 
model with exponential cutoff. Typical values for the cutoff 
power law model parameters are photon index a = 1.1 — 2.5; 
cutoff energy 5 — 18 keV; e-folding energy 11 — 26 keV. 
For any spectral model covering the range 1-10 keV, it is 
also necessary to include a gaussian component represent- 
ing iron line emission at ~ 6.4 keV, and a term to account 
for the effects of photoelectric absorption by cold gas along 
the line of sight with hydrogen column density in the range 
nu = (4—140) x 10 22 cm -2 . The source spectrum and in par- 
ticular the column density uh have previously exhibi ted sig- 
nificant va riability on time-scales as short as a day (Becker 
et al. l!)76). Measurements of spectral variation with phase 
are few; one example of pulse-phase spectroscopy was un- 
derta ken with data from the Ginga satellite from 1987 and 
1988 (Dotani et al. 198£). Only the column density and the 
iron line centre energy were allowed to vary with phase in 
the spectral fits, and no significant variation was observed. 
The Ghosh and Lamb (1979) model predicts a correla- 



tion between torque and mass transfer rate (and hence lumi- 
nosity) for accretion-driven X-ray sources. For most sources 
it is difficult to test such a relationship since the range of 
luminosities at which they are observed is limited. However 
the correlation between torque and luminosity has been con- 
firmed, at least approximately, for three transient sources us- 
ing data from the Burst and Transient Source Experiment 
(BATSE) aboard the Compton Gamma Ray Observatory 
(CGRO) and EXOSAT (Reynolds et al. 1996, Bildsten et 
al. 1997). 

The situation for persistent pulsars is, however, less 
straightforward. The BATSE data have demonstrated that 
in general the torque is in fact uncorrelated with luminosity 
in these sources. The spin-up or spin-down rate can remain 
almost constant over intervals (referred to in this paper as 
a 'constant torque state' or just 'torque state') which are 
long compared to other characteristic time-scales of the sys- 
tem, even when the luminosity varies by several orders of 
magnitude over that time. Transitions between these states 
of constant torque can be abrupt, with time-scales of < 1 
d when the two torque values have the same sign; alterna- 
tively when switching from spin-up to spin-down (or vice- 
versa) the switch generally occurs smoothly over a period of 
10 - 50 d. 

It is possible that there remains some connection be- 
tween the torque and luminosity, since at times the torque 
measured for GX 1+4 has been anticorrelated with luminos- 



ity ( Chakrabarty et al. 1997). This behaviour has not been 
observed in other pulsars. One important caveat regarding 
the BATSE measurements is that the instrument can only 
measure pulsed flux. Systematic variations in pulse profiles 
or pulse fraction could introduce significant aliasing to the 
flux data, hence masking the true relationship between bolo- 
metric flux and torque. Given that pulse profile shape and 
t orque state have shown evidence for correlation in GX 1+4 
(Greenhill, Galloway and Storey 1998) this could potentially 
be an important effect. 

In this paper we present results from spectral analysis 
of data obtained from GX 1+4 during 1996 using the Rossi 
X-ray Timing Expl orer satellite (RX TE; Giles et al. 1995). A 
companion paper dGiles et al. 199S| ) contains detailed anal- 
ysis of pulse arrival times and pulse profile changes. 



2 OBSERVATIONS 

The source was observed with RXTE between 1996 July 
19 16:47 UT and 1996 July 21 02:39 UT. Several interrup- 
tions were made during that time as a consequence of previ- 
ously scheduled monitoring of other sources. After screening 
the data to avoid periods contaminated by Earth occulta- 
tions, the passage of the satellite through the South Atlantic 
Anomaly (SAA), and periods of unstable pointing, the total 
on-source duration was 51 ks. RX TE consists of three instru- 
ments, the proportional counter array (PCA) covering the 
energy range 2-60 keV, the high-energy X-ray timing exper- 
iment (HEXTE) covering 16-250 keV, and the all-sky mon- 
itor (ASM) which spans 2-10 keV. Pointed observations are 
performed using the PCA and HEXTE instruments, while 
the ASM regularly scans the entire visible sky. 

The background-subtracted total PCA count rate for 

3 of the five proportional counter units (PCUs) comprising 
the PCA is shown in Fig. |l|a. The other two PCUs were only 
active briefly at the beginning of the observation so those 
data are not included in the analysis. The phase-averaged 



PCA count rate was initially low at 



) count s 



This 



corresponds to a flux of ~ 6 x 10 ergs" in the 2-60 keV 
energy range, using the spectral model discussed in section ^ 
and assuming a source distance of 10 kpc. Throughout this 
paper we shall use this value as the source distance unless 
otherwise specified; the actual distance is thought to be in 
the range 3-15 kpc (Chakrabarty and Roche 1997). During 



the course of the observation the count rate decreased to 
a minimum of ~ 5 counts -1 , corresponding to a flux of « 
4 x 10 35 ergs -1 before partially recovering towards the end. 
The count rates are unusually low for this source, with other 
observations giving significantly higher rates; for example 
w 320 counts -1 and 230 counts -1 (equivalent rates for 3 
PCU's) in February 1996 and January 1997 respectively. At 
times the background-subtracted countrate during interval 
2 drops significantly below zero. This is a consequence of the 
low source to background signal ratio (around 1:10) during 
this interval coupled with statistical variations in the binned 
countrate values. 

The observation is divided into three intervals on the 
basis of the mean flux (Fig. jj]a). Interval 1 covers the start 
of the observation to just before the flux minimum. Interval 
2 spans the period of minimum flux, during which time the 
flux was 30 count s -1 apart from ~ 10 s during a flare (see 
section ^|). Interval 3 covers the remaining portion of the 
observation, during which the mean flux increased steadily. 

Accompanying the changes in flux were significant vari- 
ations in pulse profile and spectral shape. Historically GX 
1+4 has shown evidence of a c orrelation between torque 
state and pulse profile shape ( |Greenhill, Galloway and 
Storey 1998j ). Throughout the period of spin-up during the 
1970s pulse profiles were typically brighter at the trailing 
edge with respec t to th e primary minimum; e.g. Doty, Hoff- 
man and Lewin (1981). Since then, measured pulse profiles 



have instead usually been leading-edge bri ght, w ith less pro- 
nounced asymmetry; e.g. Greenhill et al. (1993). During in- 



terval 1 the pulse profile was observed to be leading-edge 
bright, similar to other observations since the 1980s. Pul- 
sations all but ceased during interval 2, and in interval 3 
the shape of the profile had changed dramatically and re- 
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Figure 1. a) Background-subtracted PCA countrate for GX 1+4 
over the course of the observation, showing the division between 
intervals 1, 2 and 3 (see text). The binsize is 16 s. b) Mean PCA 
spectra between 2.2 and 40.0 keV during each interval. The data 
mode used is Standard-2, with 128 channels over the total range 
of the instrument, c) Ratio of mean spectra in intervals 2 and 3 to 
that of interval 1, labelled 2/1 and 3/1 respectively. The dotted 
line indicates the energy of the Fe line. Spectra were obtained 
from the PCA except for the point covering the highest energy 
range, which was calculated using HEXTE data. 



sembled the traili ng-edge bright pr ofiles typically observed 
during the 1970s flGiles et aL 1999| ). 

The count rate spectra taken during each interval are 
shown in Fig. |l]b. The overall spectral shape changed signifi- 
cantly over the course of the observation, with the spectrum 
becoming harder in intervals 2 and 3 compared to inter- 
val 1. The iron fluorescence line at around 6.4 keV appears 
more prominent in the second and third intervals. Iron line 
enhancement during intervals 2 and 3 is also apparent in 
the spectral ratios, Fig. [j]c. These ratios were calculated by 



subtracting the background spectrum (including a compo- 
nent to account for the emission from the galactic plane; 
see section ^|) from the source spectrum for each interval 
and dividing the resulting spectra for intervals 2 and 3 by 
that of interval 1. Because the countrate drops off steeply 
above 10 keV the spectral bins must be made correspond- 
ingly larger to achieve a reliable ratio. The datapoint in 
the highest energy band for each curve was obtained from 
HEXTE data, while the lower energy ratios are calculated 
from PCA data. The decrease in flux observed from interval 
1 to 2 and 3 becomes more pronounced at energies below 
6 keV. Above 15 keV the spectral ratios are almost constant 
with energy. 



3 THE SPECTRAL MODEL FITS 

Instrumental background from cosmic ray interactions and 
as a result of passages close to the SAA are estimated us- 
ing the PCABACKEST software, provided by the RXTE Guest 
Observer Facility (GOF) . Due to the proximity of the source 
to the galactic plane, an additional component which takes 
into account the so-called 'galactic ridge' emission must be 
included in any spectral model. The model for this com- 
ponent used in our fit s is identical to that fitted to survey 
data from this region (Valinia and Marshall 1998) with nor- 



malisations and abundance fitted to spectra taken during 
slews to and from the source during this observation. A sec- 
ondary instrumental effect which must be taken into account 
to obtain the lowest possible residuals in the model fit is a 
consequence of the Xenon L-edge absorption feature. This 
feature is modelled in our spectrum by a multiplicative edge 
model component with energy fixed at 4.83 keV. 

Candidate spectral models were tested by fitting to the 
count rate spectrum to mini mise \ 2 using t he xspec spectral 
fitting package version 10 flArnaud 199q ). In general each 
model takes the form of one (or more) continuum compo- 
nents with a gaussian component necessary to simulate the 
iron line emission, and a multiplicative component to ac- 
count for absorption by cold matter along the line of sight. 
With a lower than normal count rate for the source during 
this observation, the primary source of error is the Poisson 
statistics within each energy bin rather than any instrumen- 
tal uncertainty. Hence the systematic error in xspec was set 
to zero for the model fits. 

Thermal bremsstrahlung and power law continuum 
components both resulted in formally unacceptable values of 
X 2 for the interval 1 mean spectrum. Some improvement was 
found by fitting with various forms of power law, including a 
power law with exponential cutoff and a broken power law. 
Nevertheless, each of these models gave unacceptable values 
of reduced-x 2 : 2.93 and 1.94 respectively. An acceptable fit 
was obtained using an analytic model based on Compton- 
isation of a thermal spectrum by hot plasma close to the 
source ('compTT' in XSPEC; Titarchuk 1994), with reduced- 
X 2 for the interval 1 spectra of 1.10. The Comptonisation 
model has thus been chosen for all spectral fitting reported 
in this paper. Model parameters for spectral fits from inter- 
vals 1 and 3 are shown in Table |l[ 

Fitting this model to the interval 2 mean spectra re- 
sulted in an acceptable xt fit statistic of 0.7943, but with 
very wide confidence limits for the fit parameters. No im- 
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provement in the confidence intervals is obtained by freez- 
ing selected parameters to the mean values for the entire 
observation (e.g. To). The fit parameters are effectively un- 
constrained and as such cannot be relied upon as a measure 
of the source conditions. Additionally, the interval 2 spec- 
tra alone do not permit an unambiguous choice of spectral 
model. We cannot distinguish between cutoff powerlaw, bro- 
ken powerlaw, and Comptonisation spectral models during 
this interval on the basis of xt- A comparable fit can even 
be obtained using a model consisting of two blackbody emis- 
sion components, with fitted temperatures 1.4i;| keV and 
6.0l;5 keV (xt = 0.81; see section ^). Consequently we re- 
strict the discussion of the mean spectral fitting results to 
those from intervals 1 and 3. 

The Comptonisation model implementation in xspec 
offers a geometry switch which affects the fitted value of 
the optical depth rp. The switch can be set to model ei- 
ther disc or spherical geometries. As we will argue in section 
neither of these are strictly appropriate for the present 
situation. Consequently we performed all the fitting using 
the disc geometry, but note that fitted values of rp with the 
spherical geometry will be approximately twice as large. The 
fitted values of rp should provide an adequate comparative 
measure of the degree of Comptonisation between different 
spectra. 

The increase in line-of-sight absorption following inter- 
val 1, suggested initially by the spectral ratios (Fig. |l|c), 
is further supported by the model fits. The fitted column 
density tig more than doubles between interval 1 and 3. 
Spectral fits to each uninterrupted 'burst' of data (see Fig. 
|l]a) indicate that the increase took place smoothly over ap- 
proximately 10 hours, although significant variations in the 
fitted nn values are observed on timescales as short as 2 h. 
BeppoSAX satellite observations indicate that ng may have 
persisted at the level measured at the end of interval 3 at 
least until August 19 (Israel et al. 1998). 

The input spectral temperature kTo is consistent with 
a constant value of « 1 keV during the entire observation. 
The decrease in flux following interval 1 is associated with 
a marginally significant decrease in the fitted values of the 
scattering optical depth rp and the Comptonised component 
normalisation parameter Ac- 

The model parameters associated with the gaussian 
component, representing fluorescence from iron in the cir- 
cumstellar matter, are consistent with constant values over 
the course of the observation. The iron line centre energy is 
consistent with emission from cool matter, with no signif- 
icant change in the centre energy found between intervals. 
The iron line equivalent width (EW) increases with marginal 
significance following interval 1. 



4 PULSE-PHASE SPECTROSCOPY 

The data from interval 1 were divided into 10 equal phase 
bins, and a spectrum obtained fo r eac h phase range. The 
ephemeris is that of Giles et al. (1999), with best-fit con- 



was not possible to fit models reliably to pulse-phase spec- 
tra from interval 2 (due to the low countrate) or interval 3 
(due to its short duration). Secondly, since the evidence of 
the pulse profiles suggests conditions in the source may be 
rather different between intervals 1, 2 and 3, this seems a 
better choice than simply combining all the data. 

Each of the 10 spectra were then fitted with the model 
described in section |3| Initially, all fit parameters (barring 
those of the galactic ridge component) were left free to vary. 
Fitted values of the column density nn, input spectrum tem- 
perature kTo, and the iron line component parameters were 
all found to be consistent with those for the mean inter- 
val 1 spectrum. Confidence limits for the scattering plasma 
temperature kT were very large within some phase ranges, 
while significant variations with phase were observed only 
in the normalisation parameter Ac and scattering optical 
depth rp. To improve the confidence intervals for the latter 
three parameters a second fit was performed with all other 
parameters frozen at the fitted values for the mean interval 
1 spectrum. The resulting fit values are shown in Fig. [| 

Around the phase of primary minimum (<f) — 0.0 — 0.1) 
the fitted value of rp is significantly higher than the mean 
value, while kT is lower. The normalisation Ac is also sig- 
nificantly lower than the mean value at <f> = 0.0, but in the 
phasebin immediately following is above the confidence lim- 
its for the mean. There is little evidence for strong spectral 
variation from the mean at other phases, however we do note 
an almost monotonic decrease in the normalisation Ac from 
= 0.1 and throughout the pulse cycle. 



5 FLARE SPECTRA 

During the period of lowest flux (interval 2) a strong flare 
was observed, with the peak flux rising to almost 20 times 
the mean level during this interval (Fig. ^). The flare was 
preceeded by a modest brightening of the source which be- 
gan ~ 150 s before the flare itself and lasted ~ 60 s; a 
second pre-flare brightening began « 50 s before the main 
flare, lasting ~ 30 s. Both the flare and the pre-flare activity 
occurred within the extent of two pulse periods. From th e 
ephemeris determined for the full data set (Giles et al. 1999), 
a primary minimum would have occurred between the first 
and second pre-flares had the source been pulsing as was 
observed during intervals 1 and 3. The instantaneous flux 



stant barycentre corrected period P — 124. 36568±0. 00020 s. 
The primary minimum is defined as phase zero. Data from 
interval 1 alone were used, for two reasons. Firstly, the coun- 
trate was at its highest during that time, making the signal- 
to-noise ratio optimal compared to the other intervals. It 



during the flare peaked at « 105 count s - , compared to the 
mean rate during interval 2 of ~ 5 count s _1 . No comparable 
events occurred at other times during interval 2. During in- 
tervals 1 and 3, the significant variations between successive 
pulse profiles make it difficult to rule out flares with peaks 
having similar heights above the mean level. Certainly no 
flares with the same proportional increase in flux compared 
to the mean occurred over the course of the observation. 

The count rate during this interval was too low to obtain 
useful full-resolution spectra. Instead, low-resolution spectra 
at various times were extracted from the uninterrupted por- 
tion of data within which the flare was observed. The PHA 
ratios obtained by dividing the various spectra are shown in 
Fig. ^|b. The top panel shows the ratio of the mean spectrum 
following the flare to that preceding it (excluding the flare 
itself). Mean flux decreased by around 50 per cent following 
the flare, with no strong evidence of spectral variation. The 
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Table 1. Fit parameters for intervals 1 and 3 using the best-fitting model based on Comptonisation of soft photons 
by hot plasma. kTo is the temperature of the thermal input spectrum, kT and rp are the temperature and optical 
depth respectively of the scattering plasma, and Aq is the normalisation parameter for the Comptonised model 
component. The model also incorporates a gaussian component representing iron line emission, with Ep c the line 
centre energy, a the line width, Ap e the normalisation and EW the equivalent width. Both these components are 
attenuated by photoelectric absorption by cold matter along the line of sight, with column density njj. Data used 
is from PCA mode Standard-2. Confidence intervals are 90 per cent; fit statistic is reduced-x 2 (x 2 )- 
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second and third plots show the ratios of the pre-flare (in- 
terval A on Fig. |3p,) and flare (intervals B and C) spectra 
vs. the mean (excluding the flare). In each case there is no 
evidence of spectral variation; each ratio is consistent with 
constant PHA ratio in the range 3-20keV. The pre-flare ex- 
hibits only a modest increase in flux of around 50 per cent, 
while during the 58 s window encompassing the flare itself 
the mean flux increased by 4-5 times. 

The bottom panel shows the PHA ratio between the 
falling and rising parts of the main flare (intervals C and 
B respectively). The ratio is constant barring a broad dip 
between 6 and 12 keV. Examination of the spectra indicate 
that this dip is not due to any global change in the spectral 
shape, but rather a localised decrease in flux within that 
energy range as the flare developed. Modelling the spectra 
from the rising and falling parts using the two-temperature 
model described in section H, we find that the variation can 
be fitted best by a decrease in the temperature of the cooler 
component (x„ = 1.6) although the change in temperature 
is not statistically significant. 



6 DISCUSSION 

Re-analysed BATSE data confirm that GX 1+4 underwent 
a torque reversal from spin-down to spin-up around 1996 
August 2, approximately 10 days after the RXTE obser- 
vation ( |Giles et al. 199E ) . Since contributions to the net 
torque on the neutron star may come from both accreted 
material and magnetic stresses within the disc, it seems rea- 
sonable to suggest that changes in the magnetosphere or 
disc structure which cause the torque reversal may occur 
some time before a measurable effect is seen on the star it- 
self. We therefore suggest that the spectral and pulse profile 
changes measured during our observation are related to the 
(presently unknown) phenomenon which causes torque re- 
versals. Additional support for this connection is provided 



by the observation of dramatic pulse profile shape changes 
during the RXTE observation coupled with the previously 
noted correlation between pulse profile sha pe and torque 
state ( Greenhill, Galloway and Storey 1998[ ). Until an ob- 
servation can be made which encompasses the precise time 
a torque reversal is occurring it may be impossible to deter- 
mine more about the process. 

Comptonisation models have been used to fit spectra 
f or th is source from past observations, e.g. Staubert et al. 
(L99fJ); however it has not previously been possible to elim- 
inate all other candidate models on the basis of the \ 2 fit 
parameter. The particular model used for the spectral fit- 
t ing simulates C omptonisation in an unmagnetised plasma 
(Pitarchuk 1994), and since the available evidence points 



towards a strong magnetic field in GX 1+4 (although this 
awaits confirmation by more direct measurements such as 
a cyclotron resonance line) the model fit parameters may 
not be an accurate measure of the source conditions. It is 
likely that the principal effect of the magnetic field will be 
to make the spectral parameters dependent on the emission 
angle. Hence the model fit parameters obtained from the 
mean spectra are expected to be a reasonable approxima- 
tion of the actual values (L. Titarchuk, pers. comm.) 

Assuming that the majority of the X-ray emission orig- 
inates from a blackbody at most the size of the neutron star 
(R* » 10 km), we expect a temperature kTo ^ 0.5 keV. 
The temperature of the input spectrum kTo ~ 1 keV ob- 
tained from the model fits is consistent with this calcula- 
tion. Rough estimates of the accretion column density can 
be made based on the mass transfer rate derived from the 
luminosity, and assuming a simple column geometry. The 
accretion luminosity L acc » GM*M / R* and hence during 
interval lMfs2x 10 16 gs _1 . Assuming that the accretion 
column radius R c is some fraction / of the neutron star 
radius and the column plasma is moving at approxi- 
mately the free fall velocity ~ 0.5c, the estimated optical 
depth for Thompson scattering is ~ 0.17/ /. In general / 
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Figure 2. a) Spectral model fit parameters with pulse phase 
during interval 1 . The top panel shows the scattering plasma tem- 
perature kT; the middle panel the optical depth for scattering rp\ 
while the bottom panel shows the model component normalisa- 
tion Aq. The errorbars show the 90 per cent confidence limits. 
The fitted parameter values for the mean interval 1 spectra are 
shown by the solid lines; 90 per cent confidence intervals as dotted 
lines, b) Pulse profile from background-subtracted Event-mode 
PCA data. Two full pulse periods are shown for clarity. 



is subject to considerable uncertainties particularly given 
the over-simplistic geometry adopted here, but we estimate 
/ « 4 X 10~ 2 (e.g. Frank, King and Raine 1992) and thus 
the optical depth r w 5, close to the model fit values. 

The pulse phase spectroscopy results also show that Tp, 
kT and Ac are significantly modulated at the pulsar rota- 
tion period. Consequently we propose that the Comptoni- 
sation model provides a realistic picture of spectral forma- 
tion in this source, with scattering taking place in the ac- 
cretion column. Thus the kT parameter can be interpreted 
as the mean temperature of the accretion column plasma 
over th e region in the column where scattering takes place. 



The model normalisation parameter Ac is somewhat more 
difficult to relate to a physically measurable quantity, since 
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Figure 3. a) Event-mode PCA lightcurve (averaged over 1 s 
bins) during interval 2 showing the flare. Thick lines show the 
pre dicted times of primary minima from the ephemeris of Giles et 
al. (199£). The three intervals of interest are labelled A, B and C, 



for the pre-flarc brightening, flare rise, and flare fall respectively, 
b) The top panel shows the ratio of the mean post-flare to the 
mean pre- flare spectrum (excluding the flare itself). The second 
panel shows the ratio of the spectrum during the pre-flare increase 
(interval A in Fig. j^i) to the mean spectrum (excluding the flare). 
The third panel shows the ratio of the flare spectra (intervals B 
and C) to the mean (excluding the flare). The fourth panel shows 
the ratio of the spectrum during the flare decrease (C) to that 
during the flare increase (B). 



both the kT and rp parameters can also affect the total flux 
from the model component. 

The spectral ratios (Fig. [j]c) and the spectral fit pa- 
rameters strongly suggest that the variations in the mean 
spectra during the course of the observation are due to two 
factors. The decrease in flux which is essentially indepen- 
dent of energy is presumably a result of decreased rate of 
mass transfer to the neutron star M. This is accompanied 
by a strong increase in absorption by cold material causing 
the flux decrease below 6 keV. 

Variations in the column density uh on time-scales of 
~ 2 h have not previously been observed in this source. 
The iron line energy and the relationship between equivalent 
width and uh are consistent with the s pheric al distribution 
of matter suggested by Kotani et al. (1999), however the 



variation is much too rapid to be attributable to the negative 
feedback effect which those authors suggest regulates mass 
transfer to the neutron star in the long term. The rapid 
variation may be an indication of significant inhomogeneities 
in the circumstellar matter, or alternatively that the giant 
wind velocity is much faster than 10 km s -1 as suggested by 



infra red observations of the companion (Chakrabarty et al 



199S 



Variation in the spectral fit parameters with pulse phase 
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may provide clues to the distribution of matter in the ac- 
cretion column. The sharp dip in the pulse profiles is asso- 
ciated with a significant increase in the scattering optical 
depth tp and decrease in the Comptonisation component 
normalisation parameter Ac (Fig. ^). Such an effect may be 
observed if the accretion column is viewed almost directly 
along the magnetic axis, resulting in a much greater path 
length for photons propagating through the relatively dense 
matter of the column; essentially an 'eclipse' of the neutron 
star pole by the accretion column. Preliminary Monte Carlo 
modelling based on Comptonisation as the source of high- 
energy photons supports this as a possible mechanism (Gal- 
loway, 1999, work in progress). Accretion column eclipses 
have previously been post ulated to explain dips in pulse 
profiles from A 535+262 dCemcljic and Bulik 199c] ) and 
RX J0812.4-3114 flReig and Roche 1999D - That the plasma 
temperature kT is also low around the phase of primary 
minimum may be related to the bulk motion of the column 
plasma, since the relative velocity of the plasma in the ob- 
server's frame will depend on orientation (and hence pulse 
phase). The velocity of bulk motion is likely to be many or- 
ders of magnitude above the thermal velocity (in the plasma 
rest frame) and so may result in observable variation of this 
fit parameter with pulse phase. The asymmetry of the nor- 
malisation Ac with respect to the primary minimum fur- 
thermore points to significant asymmetry of the emission on 
the 'leading' and 'trailing' side of the pole. Such asymme- 
try may originate from a nonzero relative velocity between 
the disc and column plasma where the disc plasma becomes 
bound to the magnetic field lin es and enters the magneto- 
sphere (Wang and Welter 1981). The additional observation 



that the width of the dip decreases with increasin g energy 
may point towards a role for resonant absorbtion ( |Giles et 
al. 199^). 

The observation of a short-duration flare during the 
minimum flux period provides a further example of previ- 
ously unseen behaviour in this source. With the peak flux 
during the flare rising to almost 20 times the mean level, and 
with no other comparable events observed during interval 2, 
it is likely that the flare was due to a short-lived episode of 
enhanced accretion. The mean accretion rat e during interval 
2 can b e est imated to be ~ 2.5 x 10 15 gs" 1 (Frank, King and 
Raine ^992j ); the increased luminosity observed during the 
flare thus implies additional accretion of at least 5 x 10 17 g. In 
order to measure the instantaneous M throughout the flare 
it would be necessary to correct for the effects of anisotropic 
emission from the neutron star surface as well as changing 
observation angle with the star's rotation. Since the geome- 
try is essentially unknown, and beam patterns rather model 
dependent, this is not yet possible. 

We do however note that the delay between the start of 
the pre-flare increase ('A' in Fig. |§|) and the flare itself is w 
150 s. The relative angular velocities of the disc plasma and 
the neutron star magnetic field lines at the inner disc radius 
imply a periodicity significantly different from that resulting 
from the neutron star's rotation. From the mean interval 2 
luminosity and the estimated surface magnetic field strength 
for GX 1+4 of 3 x 10 13 G we estimate the inner disc radius 
as 2.7 x 10 7 m. A locally dense patch of plasma rotating with 
Keplerian velocity in the disc would pass close to the region 
where plasma enters the accretion column originating from 
each pole every 150 s or so. Thus it is conceivable that these 



two events represent successive passages of the same patch 
through the column uptake zone in the disc. After the second 
passage, the patch is presumably completely transmitted to 
the star and so no further flaring behaviour is seen. 

If M variations are a significant factor in the evolution 
of the flare, we might see other indications in the flare shape. 
If the polar region cools much more slowly than the flare 
timescale an asymmetric flare might be observed. Spectral 
model fits might also indicate cooling of the emission com- 
ponent originating from the pole. However, the flare appears 
almost completely symmetric, and spectral fits to the rising 
and falling parts of the flare do not exhibit cooling at any 
statistically significant level. 
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